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Introduction 

The objective of the present project is to develop a methodology of 
fracture toughness characterization of ceramics on the basis of Crack 
Diffusion Model (CDM). The latter intends to explain the scatter of 
conventional fracture parameters and the "scale effect", i.e., dependence of 
fracture parameters on the specimen dimensions. CDM recognizes the 
presence of a random field of microinhomogeneities as the source of the 
scale effect and the scatter. Statistical parameters of the random field 
characterize the material and are independent of the specimen geometry. 
Their evaluation requires a data bank of either critical loads or critical 
crack lengths or crack penetration depths. The type of the specimens is 
not restricted by CDM, but is rather dictated by convenience of testing 
procedures. 

For ceramics, such specimen types as single edge notched, bending, 
etc. used to evaluate critical quantities (load, stress intensity factor and so 
on) are undesirable. Indeed, the specimens are costly whereas each yields 
only one data point out of, say, twenty required for a statistically 
representative sample set. This motivated the choice of a compact tension 
specimen. The challenge became to achieve multiple crack initiations and 
arrests in a CT specimen made of ceramics (even single crack arrest was 
believed to be nearly impossible to obtain due to the rigidity of ceramics). 
Recently, after a year long effort, A. Calomino from NASA Lewis Research 
Center solved successfully this problem. In the meantime, the 
methodology of multiple crack initiation and arrest analysis has been 
tested on a model material (epoxy). Epoxy offers the convenience of 
specimen preparation and well detectable crack initiation and arrest 
positions. Our findings in this matter are presented below. 
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Background 


Griffith’s criterion, i.e. Gi=2y, where Gi is the energy release per unit 
crack advance (ERR) and y is the surface (fracture) energy, is one of 
the most celebrated results in studies of brittle fracture. It expresses 
the necessary condition of stability (or instability) of a cracked solid 
with respect to crack extension. Particularly, it implies that the 

(a) 

critical value of ERR at an arrest of a propagating crack, Gj , coincides 

with that at crack initiation, G^, both being equal to 2y. 

However, experimental data on brittle crack formation in epoxides 

indicate that the values of G^ may be significantly larger than those 


( a ) 

of Gj 1-5. This difference depends on loading rate, temperature, 
crack propagation rate, chemical composition and curing conditions of 
the epoxy. The origin of the difference remains to be identified. In 
this paper a systematic energy analysis of crack initiation and arrest, 
quantification of the micromechanisms of this phenomena and a 

discussion on the origin of the difference between G^ and G^ are 
presented. 


Material, specimen geometry and testing procedure 

Two plaques of thickness 2.1 mm and 4.2 mm were prepared from a 
mixture of Dow Epoxy Resin 332 (a bisphenol A diglycidyl ether), 

Dow Epoxy 732 (a polyglycol extended diglycidyl ether) and Dow 
Epoxy Hardener 20 (a diethylene-triamine curing agent) cured under 
the same conditions, (room temperature, 48 hours). The weight ratios 
of DER 332:DER732:DEH 20 were 100:10:12.5. The Tg of the plaques 
were determined to be 61±2°C using a Differential Scanning 
Calorimeter with a heating rate of 10°C/minute from 25°C to 200°C. 

Two sets of compact tension specimens (small and large with 
dimensions shown in Fig. 1) were machined from the plaques. The 
initial cracks were introduced by razor blades. The thickness of 
individual specimens is reported in Table 1 and 2. 
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Tension ramp tests were performed on an Instron 1330 servo- 
hydraulic test system at ambient conditions. All tests loadings were 
displacement controlled with a rate of 0.14 mm/min. 
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Figure 1 Specimen geometry and sizes; w = 25.6 mm and 45.0 mm 
for small and large specimens, respectively. 


Results and Observations 

Typically 4-5 crack initiations and subsequent arrests are gained 
from each specimen. Thus, just a few specimens supply a statistically 
representative set of data for analysis. A typical load-displacement 
curve is shown in Fig. 2. The load rises proportional to the increase in 
displacement to the first crack initiation point, 1*. At this moment the 
crack rapidly advances, causing a sudden drop in the recorded load 
until the crack arrests at position l a . (Here, we employ the term 
"initiation" to indicate a transition from a steady crack or a slow 
quasi-static crack growth to a fast dynamic crack advance). As the 
displacement continues this process is subsequently repeated for a 
number of steps with crack initiations depicted at positions 2^,3*, etc. 
and corresponding arrests at 2 a , 3 a , etc.. This phenomena is 
conventionally known as "stick-slip". 



Load 
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The tangents of the load-displacement curves at the points l a , 2 a , 
depicted by the dotted lines in Fig. 2, go through the origin, 
indicating a perfectly elastic response of the system immediately 
after crack arrest. On finer observation of the load-displacement 
curve on increasing displacement after crack arrest, shown in the 
enlargement of step 2 a -3‘, the curve deviates from linearity (i.e. an 
increase in compliance) and consists of smaller sharp drops of load 
resembling the larger steps. The increase in compliance indicates a 
sub-critical crack growth and the saw-tooth nature of this curve, 
shown in the enlargement, that this process procedes by discrete 
events. These suggestions are well supported by the morphology of 
the fracture surfaces. 

Figure 3a shows a sketch of the fracture surface of specimen 2 
observed by optical means. The dotted and solid lines indicate crack 
arrests and initiations, respectively. Triangles, or "deltoids", form 
between the arrest and initiation lines during sub-critical crack 
growth. The fact that the deltoids do not appear before the arrest 
line was established by rapid fracture of a specimen cooled in liquid 
nitrogen after the crack arrest. Figure 3b is an optical micrograph of 
the deltoid formations between the second arrest and third initiation 
events in specimen 2. These deltoid formations on the fractured 
surfaces are common for all specimens examined here and also can 
be seen in other studies 3 * 4 . The deltoids were further examined by 
scanning electron microscopy, SEM, after platinum coating the same 
specimens previously inspected by optical means. Figure 3c shows 
an SEM micrograph at 1000X of the region designated by the small 
box "B" in Fig. 3b. The deltoids appear on the scale of microns in 
height and 100 microns or larger in width and are complimentary on 
the two crack faces. The complimentarity, which has also been 
observed elsewhere 4 , excludes any inelastic deformation from 
consideration. In other studies with epoxies prepared under various 
conditions some localized plastic shear deformation in the vicinity of 
the crack tip has been reported 1 * 2 * 5 . However, Kinloch et. also 

reported the complimentarity of the fracture surface which should 
not be observed if plastic deformation occurs. It should also be noted 
that according to Phillips et al 4 the results of previous authors could 
not be reproduced (see p. 322 in ref. 4). The sub-micron scale texture 
of the fracture surface of the deltoid is similar to that of the crack 
jump after initiation and shows no indication of ductility. 


Arrest line 
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DIRECTION OF CRACK 



[ b ] [ C ] 

Figure 3a A sketch of the fracture surface of specimen 2. 

3b An optical micrograph of the deltoid formations between 
the second arrest and third initiation events in specimen 2 
3c A SEM micrograph at 1000X of the region designated by 
the small box "B" in Fig. 3a. 
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Figures 4a and b show representative optical micrographs of the sub- 
critical crack growth zones of a small and large specimen, 
respectively, at the same magnification. The deltoid formations 
appear as a cascade of elementary deltoids with their bases 
extending in the direction of crack growth. The size of an elementary 
deltoid clearly appears to be invariant with respect to specimen 
thickness and also with the steps l a -2‘, 2 a -3‘, etc. Thus, it appears 
that the discontinuous nature of the sub-critical crack advance 
recorded in the load-displacement curves is related to the deltoids 
observed on the fracture surface. 



a 


b 


Figure 4 Optical micrographs of the cascades of deltoids of a small 
[a] and large specimen [b], at the same magnification. 
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Discussion 


The energy release rates at initiation, G/ , and arrest, G* , as well as 
the total energy released on sub-critical growth, can be readily 
evaluated based on the elastic solution (Appendix 1). The values of 

the energy release rates G^ and G (a) for specimens 1 and 2 are given 
in Figs. 5 and 6. Also shown are the corresponding fracture surfaces 
between the crack arrest and initiation lines. Typically, the values 


for G (a) for all the specimens are fairly constant and significantly 

smaller than those for G^\ The values of G^ change randomly with 
crack length. Analogous findings have been reported by Phillips et 
al 4 . The data necessary for ERR evaluation for the specimens studied 
here are summarized in Tables 1 and 2. 

It is our observation that the ERR at initiation is strongly correlated 
with the number N of deltoids in the last row of deltoid cascades, 
before initiation of the crack jump, as shown in Fig. 7. The 

intersection of a linear regression of G^ ^ as a function of N, in the 


limit of N approaching zero, falls within the range of 



(0.16 KJ/m 2 ) 


in agreement with Griffith's criterion. The tendency of G^ to 
( a ) 

approach Gj with increase in rate of loading as well as decrease of 
temperature have also been observed 1 * 2 ’ 4 ' 6 . 


It is expected that the larger the difference Gj' - G* i.e. the larger 
the excess of energy release rate at initiation over that at arrest, the 
larger the subsequent crack jump. An estimation of the crack jump, 
Al, can be readily obtained within the framework of linear-elastic 
fracture mechanics. It should be noted that the crack jump after the 
initiation 1M a , 2^-2 a , etc., (Fig 2) is a faster process than the 
response of the actuator of the loading system. Thus the 
displacement, 8i, at crack initiation remains practically constant 
during the crack jump from l[ to l[+Al. After the jump, the crack 


arrest Gr ; and initiation G 



Energy release rate at crack arrest and initiation vs. 
normalized crack length, // w, and the micrographs of the 
fracture surface between the corresponding arrest and 
initiation lines of specimen 
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Figure 6 Energy release rate at crack arrest and initiation vs. 

normalized crack length, // w, and the micrographs of the 
fracture surface between the corresponding arrest and 
initiation lines of specimen 2. 
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Table 1 Small specimen 


































































































































































arrest takes place when the current energy release rate is reduced to 

(a) 

the value at arrest, This can be expressed as follows: 


Gi(; i + A(,8i)=G < 1 a) 


( 1 ) 


The LHS of the equation 1 can be expanded into a Taylor series 


9G, 


G 1 (/ i +A/,8 i )=G 1 (/ i 8 i )+— Ij.g. A / + 

9/ v 1 


i roi 

r i / i.8i 


a r + ... 


9/ 


( 2 ) 


Neglecting the omitted higher order terms in the equation 2 and 
taking into consideration the equation 1 and the condition of crack 

initiation, Gi(/i ,8i)=G^ , one obtains the following quadratic equation 
for A/: 


1 a 2 Gj 

2 9/ 2 


Mi 


A/ 2 + 


9Gj 

9/ 


Mi 


M 


+ G ( 1 i) -G ( 1 a) = 0 


(3) 


The solution of the equation 3 normalized by the specimen width, W, 
is shown in Fig. 8 as a function of the normalized crack length, // W. 
The dotted lines represent the upper and lower limits of the crack 
advance corresponding to the minimal and maximal values of 


,(i) 


40 


(G! - Gj )/ Gj . The solid line in Fig. 8 discribes the solution of the 

equation 3 for the average value of (G^ - G^)/ G^ . Figure 8 shows 
that the experimental data points are all within the band expected 
from the elastic solution. The scatter of the data in the band 
apparently reflects a mechanism controlled by chance. 


Utilizing the same elastic solution we estimate the total potential 
energy release, An, due to sub-critical crack advance (Appendix 2). 
A new surface area created during this process is A/ t (t is the 
specimen thickness) and the part of this area, N A D , is occupied by a 
cascade of deltoids (A D is the elementary deltoid area and N is the 
total number of deltoids in cascade). Thus an excess of An over 
2y 0 A/ t can be attributed to the formation of deltoids: . 
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Figure 7 Energy release rate at initiation vs. No. of deltoids in the 
last row of the corresponding cascade. 
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Figure 8 The normalized crack advance A//w between the 

initiation and the following arrest line vs. // w. The solid 
and dotted lines represent average minimal and maximal 
values, repectively. The data points are observed values. 
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(4) 


1 5 


All - 2y 0 A/t = 2(yD - y 0 )N-A D 


Here 2y 0 = Gr, a) stands for the specific fracture energy (SFE) for a 
smooth fracture surface and 2y D reflects the SFE for a fracture 
surface of deltoids. Figure 9 is a graphical representation of the 
equation 4 normalized by the crack advance area, A/t; specifically it 
shows the potential energy release per unit area of sub-critical crack 
advance vs. the normalized deltoids area, (N-A D )/(A/ t). 

The dashed line represents a linear regression of all the data points 
with a fixed intersection at the average value of the ERR at crack 
arrest. 

(a) 

The intersection corresponds to Gj and the slope reflects the excess 
of 2y£ over 2y 0 . The average of the deltoid SFE, <2y D >, is 0.261 KJ/m 2 
and the standard deviation is 0.031 KJ/m 2 . 

This excess of energy may be associated with kinetic, heat and/or 
other forms of energy. 
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Figure 9 The potential energy release per unit area of sub-critical 
crack advance vs. the normalized deltoids area, N A D /A/ t. 




From the observations, the creation of the deltoids appears to be 

controlled by chance, which results in a large scatter of G^. The 

difference between the fluctuating G^ and relatively constant (j** i.e. 
AG i is evaluated for the small and large specimens. The results reveal 
a strong scale effect, depicted in Fig. 10 which shows the 
distributions of AGj fitted by gamma distribution. 

The scale effect resembles that expected for strength of a material 
controlled by the largest defects. This resemblance, together with the 
characteristic appearance of the apex of the deltoids, depicted in Fig. 
4, are strongly suggestive of an underlying defect population as a 
primary source for the observed phenomena. Indeed, the deltoids 
are analogous to those created on a river surface downstream from a 
small obstacle. 

However, in spite of our expectation based on this analogy and the 
observed scale effect, there are no "defects" at the apex of the 
deltoids detectable down to the 0.1 micron level. Moreover, the 
arrest and initiation lines do not correlate with apparent , large-scale 
defects such as bubbles. This suggests . that the formation of deltoids 
during the slow process of sub-critical crack growth is an intrinsic 
feature of the fracture process itself. It may be triggered by 
inhomogeneity of the epoxy. The inhomogeneity could be related for 
instance to a fluctuation in the cross-link density. A number of 
studies on epoxies have reported nodular structures or chemical 
crosslink density fluctuations which could arise out of inadequate 
mixing of the highly viscous components 5 * 7 ' 9 . The relationship among 
the scale of structural inhomogeneity, the loading history and the 
micromechanisms of failure would thus be important to uncover. 



f ( A G 



Figure 10 Distribution of AG i =Cj \ * - fitted by Gamma 
distribution for small and large specimens. 
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Conclusion 

1. Four-five crack initiations and arrests are gained from each 
specimen. Thus following the experimental procedure described in 
this study, three specimens supply a statistically representative 
data. 

2. The deltoids are formed on the fracture surface during the slow 
sub-critical crack growth, prior to the smooth mirror-like surface 
of the fast crack. The two faces of the fracture surface are 
complimentary which excludes any inelastic defomation. 

3. The values of ERR at arrest are fairly constant and significantly 
smaller than those at initiation. 

4. The values of ERR at initiation display random scatter and strong 
correlation with the number of deltoids in the last row of deltoid 
cascades. 

5. There is a well pronounced scale effect in ERR at initiation. It 
resembles that in conventional strength of a material controlled by 
the largest defects, though no defects are found at deltoid apexes. 


Acknowledgements 

The financial support of this work by the NASA Lewis Research 
Center (Grant No. NAG3-1034) and partial support by the Office of 
Naval Research (Grant No. ONR N00014 J 1306) are gratefully 
acknowledged 



Appendix 1 


The Stress Intensity Factor (SIF) for CT specimens shown in Fig. 1 is 
evaluated as follows(ASTM 399-81 10 ) : 


S.I.F. K r 


g E ' (2+%)(0. 866+4.64%- 1 3.32q + 14.7 2q -5.601;) 
V©Vw 


where 1 1 


V(4) = 


1.277-t- 15.45^+1.33sin(2II^) 

(i-4) 2 


d -«2 


for 0.2<^<0.8 


A-l.l 


A-1.2 


and ^=//w The effective Young's modulus can be directly evaluated 
from the load-displacement curves such as shown in Fig. 2. Then the 
energy release rate Gj = Kj 2 / E' can be readily evaluated as a 
function of the crack length / and the load point displacement (LPD), 

5. The dependency of Gj on the normalized crack length £=//w for 
different LPD, 8, is shown in Fig. 11. 



G l/E' W ( x Mr 


20 



/ / w 


Figure 11 A family of normalized ERR vs. normalized crack length, 
// w, parameterized by the normalized load-point 
displacement, 8/w. 




21 


Appendix 2 


The evaluation of total potential energy change due to sub-critical 
crack growth can be done as following: 


Since the potential energy n(P, /) or I1(P, 8) has two 
variables, the differential of IT is determined by: 


jrT ari an 

dll = — -d/ + — — dP 

di ap 


independent 


If we approximate 5=8(P) from the experimental results, then the 
elastic equilibrial crack / (P, 5) becomes a function of P only. 
Therefore 


d/=^dp 

dP 


and consequently A-2.1 becomes 

JTT fan di an\ „ 
dn ” ( a/ dP + ap J dP 


where 


A-2.1 


A-2.2 


an 

a/ 



an i 

aP 2° ; G\ = Gj x (specimen thickness) A-2.3 


If one approximates 8=8(P) by a parabolic function 
5 = 5 k +A(P-PD-B(P-PD 2 


A-2.4 


where A and B are constants and can be determined with such 
conditions: 


d£ 

| j^p, = C(/ k a ) ; C(/ k ) stands for the compliance 
6 k+ i=5 k +A(Pi +1 -P|) - B(Pi +r P|) 2 
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After the calculation 


A = C«), B=— ® 


J k+1 


5 k 


<p k+ i-pj) 2 


A-2.5 


From A-2.2 and A-2.3, the total potential energy change 


An = J “ [ - Gj{/( P),P)}^- ~y5(P) ]dP 


The first term of A-2.6 


A-2.6 


M d Gj {/(P),P) 


dP 


pL ^ pi 

j . ["G 1 {/(P),P}^]dP = -G 1 {/(P),P}/(P)| p^ + J. /(P) 

= J^Tk+l + ^lk)^k+l _ ^k) 


The second term of A-2.6 

_ (Pk + 1 - Pk) r {C(4 +1 ) + C(4)KP u t+pj) 1 

" 2 L 4 


-dP 


A-2.7 


+ T {C(4)(Pj c+1 -P^)-(6 k+1 -8 k )} A-2.8 


] 


Thus Eq. (A-2.6) = Eq. (A-2.7) + Eq. (A-2.8) 

and An is expressed through observed quantities as 

Pk . pjckl . 5 k . S k+ 1 . Gfk • G ( ,i, , C(© , C(4 +1 ) , 4 , 4 +1 
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For the values of parameters under consideration the first term is 
about 95% of An. Therefore the total potential energy change can be 
approximated as follows: 

An = jWLi+CaXCi-'i) 
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